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Abstract. A new analytic solution for 2nd-order Fermi acceleration is presented. In partic- 
ular, we consider time- dependent rates for stochastic acceleration, diffusive and convective 
escape as well as adiabatic losses. The power law index q of the turbulence spectrum is 
unconstrained and can therefore account for Kolmogorov {q = 5/3) and Kraichnan (q = 3/2) 
turbulence, Bohm diffusion (q = 1) as well as the hard-sphere approximation (q = 2). This 
considerably improves beyond solutions known to date and will prove a useful tool for more 
realistic modelling of 2nd-order Fermi acceleration in a variety of astrophysical environments. 
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1 Introduction 

Stochastic acceleration of relativistic particles by plasma wave turbulence ~ a 2nd-order Fermi 
process [1] - is dominating the production of non-thermal particle distributions in a variety 
of astrophysical environments. Radio galaxies [2-7], clusters of galaxies [8-11], gamma-ray 
bursts [12, 13], extra-galactic large-scale jets [14, 15], blazars [16, 17], solar flares [18-20], 
the interstellar medium [21, 22], the galactic centre [23, 24], supernova remnants [25-27] and 
even the recently discovered "Fermi bubbles" [28] have been suggested as sites of stochastic 
acceleration. 

The dynamics of the phase-space density /(p, x, t) of relativistic particles interacting 
with a turbulent magnetised plasma is governed by a Fokker-Planck equation. If the time 
for pitch-angle scattering is much smaller than the timescales of interest, i.e. acceleration, 
escape and loss times, it suffices to consider the isotropic part of the phase-space density, 
f{p, x, t) where p = \/p^. Here we consider only relativistic energies such that energy E 
and momentum p are essentially the same, E = pc. Furthermore, we constrain ourselves to 
a distribution function f{p,t) independent of position, e.g. the spatial average if the rates 
only change slowly over the acceleration region. In the special case of transport coefficients 
constant in time (see discussion in Sec. 2), this constraint can be relaxed as the spatial and 
the momentum parts of the problem decouple, see e.g. Ch. 14 of Ref. [29]. 

Stochastic acceleration is a biased diffusion in momentum space and results in a broad- 
ening and systematic shift of the injection spectrum to higher momenta. In quasi-linear 
theory, particles interact resonantly with turbulent plasma waves of a range of wave-lengths 
similar to the particles gyroradius r^, and the resulting momentum diffusion coefficient Dpp 
depends on the spectrum of turbulence. In particular, if we consider resonant interaction 
with a power law spectrum yV{k) cx k~'^ of MHD waves of velocity = /3ac, the diffusion 
coefficient takes the form [30-32] 



pp ~ 



'9 ^2 
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where ( = {6B)'^/B'^ is the energy in turbulence {SB)"^ = J^^ dk'W{k'), compared to the 
energy of the background magnetic field B and A2 = 271 /ki is the longest wave-length of 
the MHD modes. The momentum dependence of the diffusion coefficient reflects the power 
law behaviour of the turbulence spectrum, Dpp (x p'^. In the Kolmogorov and Kraichnan 
phenomenologies the spectral index is q = 5/3 and 3/2, respectively [33]. A well known 
and particularly straightforward approximation is the so-called "hard-sphere" limit in which 
Dpp (X p^. If one however assumes the scattering mean- free path to be equal to the gyroradius 
one finds Dpp oc p (Bohm limit). 

In eq. 1 . 1 we have omitted a numerical factor which depends on q as well as on the mag- 
netic and cross helicity, e.g. in the case of slab-Alfven turbulence [34]. For 1.5 ^ q ^ 2.5 and 
the simplest case of right-/left-handed waves propagating parallel and antiparallel with the 
same power, this factor is 0(1). If Alfven waves however only propagate into one direction, 
the momentum diffusion coefficient vanishes. Furthermore, also allowing for obliquely propa- 
gating modes one needs to consider that Alfvenic turbulence is inherently anisotropic [35, 36] 
which reduces the efficiency of second-order acceleration. In contrast, fast-mode waves or 
turbulence at super- Alfvenic scales (which cascades hydrodynamically) are isotropic [36] and 
probably dominate 2nd-order Fermi acceleration [37, 38]. 

Despite the relevance of stochastic acceleration for a variety of astrophysical environ- 
ments, however, only solutions for a somewhat limited range of q and assumptions about the 
acceleration, escape and loss rates have been presented. The first solution of the Fokker- 
Planck equation including Bremsstrahlung losses but only a constant escape rate was given 
for the hard-sphere limit {q = 2) [39] using Mellin transforms with respect to particle en- 
ergy. Still for g = 2, this was later extended [40] to include time-dependent escape but 
no Bremsstrahlung. A comprehensive review of known time-dependent solutions (and a 
discussion of the importance of boundary conditions and its relevance for the steady state 
spectrum) was presented in Ref. [41], although only for time-independent rates and limited 
values of q. Finally, the transport equation has been solved for arbitrary q [29] and exten- 
sively discussed [42, 43]), however, again only for time-independent rates. We conclude that 
for time- dependent rates, a solution only exists for q = 2 whereas for general q the loss and 
gain rates must be assumed to be time-independent. 

For certain classes of environments, however, this assumption is difficult to justify. 
For example for very young (< 100 yr) supernova remnants, the acceleration, escape and 
adiabatic loss rates are expected to vary on timescales less than or equal to the acceleration 
or loss times [26, 44]. Another example are blazar jets with variability on timescales down to 
minutes [45]. What is therefore needed is the time-dependent solution of the Fokker-Planck 
equation for arbitrary q and allowing for time- dependent energy gain and loss rates. In 
this paper, we present a time-dependent solution of the Fokker-Planck equation for general 
q considering time- dependent stochastic acceleration, escape and adiabatic losses or gains. 
We employ a combination of integral transforms to reduce the transport equation to the 
heat equation in (7 + l)-dimensional spherical coordinates (where 7 is a function of q). Our 
result improves beyond solutions known to date and constitutes an important contribution 
to particle transport theory. 

In Section 2 of this paper we derive the Green's function of the Fokker-Planck equation 
with time-dependent stochastic acceleration, diffusive or convective escape and adiabatic 
loss/gain terms for arbitrary q, using a combination of integral transforms. We apply this 
newly found solution to four specific (toy) models of time-dependencies in Section 3. In 
particular, assuming all rates to be constant the solution of Ref. [42] is recovered which 



- 2 - 



constitutes a non-trivial test of our calculation. We conclude in Section 4 with some remarks 
on boundary conditions and the existence of the steady state solution. 



2 The Green's function 



We start from the transport equation for the isotropic and spatially averaged phase space 
density f{p,t) in flux conservation form [32, 46], 

m..t)_ 1 _n (,,,,?VM + ^(,,,);(,,,,^V4^ + ^. (2.1) 



dt p"^ dp \ \ ^ dp ' ' / / t{p, t) 4:TTp 

The terms on the right hand side describe biased diffusion in momentum space with a diffusion 
coefficient Dpp{p,t), additional energy gain/loss processes with a rate A{p,t), global escape 
with a rate 1/t(p, t) and injection with a rate S{p,t) / {Airp^). 

For adiabatic losses/gains the rate is proportional to p, so we define a{t) by 

A{p,t) = mc(^)a{t), (2.2) 

where m is the mass of the particle. We note that this form can in principle also account for 
bremsstrahlung losses and gains by Ist-order Fermi acceleration at shocks. However, with 
this particular momentum dependence, cooling by synchrotron radiation or inverse Compton 
scattering cannot be accounted for because the momentum dependence of the loss rate A{p^ t) 
is more complex then the usually assumed p^ which is only valid in the Thomson regime. 
While in certain limits, e.g. in the steady-state case (see e.g. [29]), analytical solutions 
might be possible, the fully general case is only amenable to numerical approaches. If the 
diffusion in momentum space is due to resonant interactions with MHD waves, the momentum 
dependence of the diffusion coefficient Dpp reflects the spectrum of the turbulence cascade. 
In particular, assuming the spectral energy density W{k) to be oc A:"*?, we have Dpp oc p'^, 
and we define the acceleration rate k{t) by 

Dpp{p,t) = k{t){mcf , (2.3) 

\mc/ 

where q = 1 for Bohm diffusion, q = 3/2 for Kraichnan turbulence, g = 5/3 for Kolmogorov 
turbulence and q = 2 in the hard-sphere approximation. 



2.1 Diffusive escape 

Assuming that the interactions with the same turbulent MHD waves dominate spatial dif- 
fusion and therefore the diffusive escape from the acceleration region, fixes the energy de- 
pendence of the escape rate. In particular, the time for diffusive escape from a region of 
spatial extent L is T -^^/ ^xx where th.6 spciticil diffusion coefficient Dxx 

is related to the 

momentum diffusion coefficient Dpp by Dxx Dpp = iv\p'^ with vx the Alfven velocity. Here, 
^ is a factor that depends on q as well as the magnetic and cross helicity of the magnetic 
turbulence. In this case, r oc p'^"^ and we define T^{t) through the relation 

r{p.t)=TS)(—Y'^ . (2.4) 
\mcJ 
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We are now looking for the Green's function to the transport equation 2.1, that is f{p, t) 
for mono-energetic, impulsive injection S{p,t) = 6{p — po)6{t — to)/{4TTp^). Introducing the 
dimensionless momentum variable x = p/{mc), the transport equation reads, 



dl 
dt 



+ 



+ 3 a{t)f + {a(t) - (2 + q)k{t)x'^-^x) x^ - k{t) 



dx"^ 



f 



6{x - xo)5{t - to) 
{mc)^4iTXQ 



(2.5) 



with xq = po/ (mc) . 

We make the substitutions, 

p{x,t) = 2x(2-«)/V£'(i)^(*) where g{t) 



exp 



{2-q) [ dt'a{t' 

Jto 



r] = ip{t) and f = f exp 
where we choose 



ya{t)- dt'X{t')a{t') -3 dt'a{t') 



to 



to 



da 
dt 



{2-qfa{tfk{t)g{t) 



il){t) = exp 



TA{t)g{t) ' 

2{2-qf f dt'a{t')k{t')g{t') 

Jtn 



^{t) = {2-qf / dt'k{t')g{t'm'). 



to 



Eq. 2.5 then transforms to 
df d^f , 1 df 



2 ' 7t^+*^^P 



dr] dp^ p dp 



-xa+ I dt'Aa + 3 / dt'a{t') 

Jto Jto 



(2.6) 
(2.7) 

(2.8) 
(2.9) 
(2.10) 



{^{t))-H{x-xo)5{t-to), (2.11) 



that is the heat equation with spherical symmetry in (7+l)-dimensional spherical coordinates 
where 7 = {4: + q)/{2 — q). Equation 2.8 is a special case of the Riccati equation and solutions 
a{t) for explicit k(t),g(t) and Td(t) are known and have been compiled in Refs. [47-49]. 
The bounded Green's function, i.e. the solution to 



that remains finite for all p and r/ > 0, is 



fiP, Po,V,Vo) 



Po 



2{ri - rjo] 



exp 



P" + Po 
4{r] - r/o) 



J 7-1 
2 



PPo 



2(r? - Vo] 



(l-7)/2 



(2.12) 



(2.13) 



with /(^_i)/2 the modified Bessel function of the first kind. Resubstituting for p = p{x,t) 
and T] = r]{t) one finds for g / 2, 



/ 



1 



2-q 



exp 



X exp 



(mc)^47rxQ 2:0 

x'^-'^gip + xI~''tPo 



dt'a{t' 



to 



exp 



2-g 2-q 
X ^ga — Xq Oq 



2-q 



Il + q 

2-q 



2(xxo)'2'' ^/g^ 



-3/2 



(2.14) 
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which, together with eqs. 2.8 - 2.10, constitutes the main result of this paper. 

For the hard-sphere approximation, one needs to carefully take the limit g — t- 2. After 
some tedious algebra including a number of non-trivial cancellations, one arrives at the known 
result [26, 40], 

2" 



/ 



1 



1 1 exp 



It 



to Td 



exp 



In X — In xo — f^^ dt'a — 3 f^^ dt'k 



(2.15) 



We stress that this result has been derived in a way completely independent from the one in 
Refs. [26, 40] and therefore constitutes a valuable test of our calculation. 

2.2 Convective escape 

If we assume an energy-independent escape time as is for example the case if convection out 
of the acceleration zone is dominating over diffusive escape, 

Tesc{p,t) = T^{t), (2.16) 

the Green's function takes a form similar to eq. 2.14, 



/ 



1 



{mcj^inxQ xq 



exp 



dt' 



X exp 



x^-'ig + xl'" 



((?-2)2/;dt'% 



to 



Il + q 

2-q 



exp 



dt'a{t'] 



to 



2-q 

(XXq) 2 ^ 

iq-2)^Jldt'kg 



2-g 



-3/2 



(2.17) 



_{q- 2)^11 dt' kg _ 

For q = 2, the result is again eq. 2.15 since in the hard-sphere approximation the escape 
time for diffusive and convective escape are both energy-independent, Tc ~ tj. 



3 Examples 

3.1 Constant acceleration and escape rates, but no adiabatic losses 

We assume that the acceleration and escape rates are constant, k{t) = ko, rd(t) = t^q, and 
that there are no adiabatic losses, a{t) = 0. The Riccati equation 2.8 is solved by 

1 1 



a 



and we find for ijj{t) and (p{t) 

m 



exp 




-2(2-(/)W^(t-to) 



(3.1) 



(3.2) 



(3.3) 



For / 2, the particle density n(x, xq, to) = 47rx^/(x, xq, t, to) for impulsive injection reads 



n(x,xo,t,to) 



2-q jl^ 2(xxo)(2"5^/^^/0 
xo \l xo {2 - q)^/kr^^{l - i)) 

^ r 4(xxo)(^~^)/^VV^ 
5^ {2 - q)y/ky^{l - 7p) 



exp 



(x 



2-9 



{2-q)Vk^il-^P) 



(3.4) 
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and for 5 = 2 this reduces to 



n{x,xo,t,to) 



1 1 exp 



t-to 



exp 



(Inx — InxQ— 3kQ{t — to)y 
4/co(i-io) 



(3.5) 



Equations 3.4 and 3.5 are identical to eqs. 46 and 49 of Ref. [42] (setting their a = 0), 
respectively (see also [29]). Their result was however derived assuming constant acceleration 
rate and escape time. Reproducing this result as a special case of our more general solution 
therefore constitutes a non-trivial test of our calculation. 

In the left column of Fig. 1 (compare with Fig. 2 of Ref. [42]) we show the particle 
spectrum for impulsive injection, n{x, l,t,0), for Kraichnan turbulence, Kolmogorov turbu- 
lence and the hard-sphere approximation. We fixed Tdo = 1 and all timescales (rates) are in 
units of l/ko (/cq). Diffusion and advection in momentum space lead to a broadening of the 
spectrum and monotonous increase of the mean energy with time. Since in the hard-sphere 
approximation the acceleration timescale p'^/Dpp is the same for all momenta, the spectrum 
(in the logarithmic momentum variable, logx) is even with respect to the mean logarithmic 
energy, (logxo + Sko^t — to))- For Kraichnan and Kolmogorov turbulence the behaviour is 
qualitatively different as the acceleration time is now increasing with energy. At higher en- 
ergies, the spectrum becomes gradually softer and asymmetric. For a fixed escape time, the 
spectrum rolls over at much smaller energies than in the hard-sphere case. 

In the right column of Fig. 1 we show the spectra for steady injection, 



n{x,xo,t,to) = dt' f{x,xo,t,t') , 



(3.6) 



to 



including the steady state spectrum, i.e. the limit t — )• oo. In general, it is not clear whether 
this integral in fact converges for t ^ oo but for the case of constant rates, injection and escape 
exactly balance each other. We note that for the hard-sphere approximation, acceleration 
and escape rates have the same momentum behaviour and consequently there is no preferred 
momentum scale. This leads to a power law steady state spectrum whereas for Kraichnan 
and Kolmogorov phenomenology the spectrum exhibits a long exponential roll-over with a 
characteristic momentum defined by the equality of acceleration and escape times. 

3.2 Constant acceleration rate 

We use the relation between the space and momentum diffusion coefficients, D^xDpp = S.v'jjP', 
to express the escape time T{t,p) = Td{t)x'^^'^ in terms of the acceleration rate k{t), 

Tdit) k{t)L^ k{t)pm \mj ' ^ ■ ^ 

where L{t) is the size of the acceleration region and the Alfven velocity VA{t) = B[t) / yfp^ 
with B the background magnetic field and pm the thermal gas mass density. Specifying the 
adiabatic loss/gain rate, a[t), to the case of an expanding/contracting flux tube [26], 



a{t) = 
with Lq = L{to 



1 f din L{t) 
3 



dt 
and Bq 



din B{t) 
dt 



g = exp 



dt'a{t') 



to 



Lq Bjt) 



B{tQ), we can write 

1 ^ 1 e 

Tdit) k{t) 

Pm 



Bo 
Lo 



9' it) 



(3. 



(3.9) 
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Figure 1. Particle spectrum n(a;,l,<:, 0) for impulsive (left panel) and steady injection (right panel) 
for Kraichnan turbulence (first line), Kolmogorov turbulence (second line) and the hard-sphere ap- 
proximation (last line) and assuming a = and t^q = 1. The solid lines are for fixed times 
i = 0.01, 0.03, 0.1, 0.3, 1, 3 and the dashed lines denote the steady state spectrum. All timescales 
(rates) are in units of 1/fco (fcg). 



The Riccati equation 2.8 now reads, 



— = a iq-2) kg — - 

dt pm \LoJ k 



(3.10) 




Figure 2. Particle spectrum n{x, 1, t, 0) for impulsive (left panel) and steady injection (right panel) for 
Kraichnan turbulence and assuming a ~ aocxp [—Xt] with oq = —1, A = 0.5 and {^B^/L^pm) = 0.1. 
The solid lines are for fixed times t = 0.01,0.03,0.1,0.3,1,3 and the dashed line denotes the steady 
state spectrum. All timescales (rates) are in units of 1/fco (^o)- 



For a constant acceleration rate, k{t) = /cq, this is solved by 




«W = «o = ^a/ — f^^' (3.11) 



leading to 



V'(i) = exp 



2{2-qfaoko [ dt'g{t') 
Jto 



(3.12) 



V{t) = :^{i^-l)- (3.13) 
2ao 

In Fig. 2 we show the particle spectrum n{x, l,t,0) for impulsive and steady injection 
and Kraichnan turbulence. We have chosen L{t) and B{t) such that 

a(i) = ao exp [-At] , (3.14) 

with ao = —1, A = 0.5 and {^Bq/ L^pm) = 0.1. Again, all times (rates) are understood to 
be in units of I/Zcq {ko). For early times {t <^ 1), the acceleration and adiabatic loss rate 
are similar, \a{t)\ ~ ko and escape can be neglected. For intermediate times {t ~ 1), the loss 
rate starts declining which leads to a rather hard spectrum. Finally, for late times {t ^ 1), 
escape becomes important such that the spectrum does not extend to any higher energies. 
The late onset of escape also reflects in very hard spectra for steady injection and a steady 
state spectrum extending over several orders of magnitude. 

3.3 Constant adiabatic loss rate and exponentially decreasing acceleration rate 

We assume that the adiabatic loss rate a{t) < is constant, a{t) = ag. This is a fair assump- 
tion for environments with a blast wave where the shock radius and speed are r{t) oc t^/^ and 
v{t) oc such that the adiabatic loss rate a{t) oc (r'^v) oc t^^^ has a weak time dependence. 
We find for g{t), 

g = go exp [At] with A = — (2 — q)ao > and go = — Ato • (3.15) 
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Figure 3. Particle spectrum n{x,l,t,0) for impulsive (left panel) and steady injection (right panel) 
for Kraichnan turbulence. We assumed a = oq = const, and fc = fco exp [—Xt] with ao = —0.6 and 
{S,Bl/Llp,n) = 0.2. The solid hnes are for fixed times t = 0.01, 0.03, 0.1, 0.3, 1, 3. All timescales (rates) 
are in units of 1/fco (fco). 



We assume further that k{t) is of the form /cq exp [—i^t] and that k = X. (In fact, this 
can be easily extended to k / A but for clarity we here constrain ourselves to the case k = X.) 
The Riccati equation 3.10 now reads 

^ = a^q - 2fko9o " - f V f exp [2Xt] . (3.16) 
and substituting y = exp [2 At], w = —1/a leads to the standard form 



A'w'+B- with ^- and B^{q-2r^. (3.17) 

dy y Pm \LoJ ko 2A 2A 

The solution for w{y) is [49] 

w = -\-^ where u = ^Ji (iVAB^) , (3.18) 



2 



Au dy 

with Ji the Bessel function of the first kind. For a{t) and tp{t) we find 



rj Ji (2Vab^) / Ji (2VABVy^, 

«(^) = V-^^/y ) ^ ( V(0 = \ ^ / I , (3.19) 

V ^ Jo (2VAB^) \ Ji (2^a^^ 



where yo = y{to). f{t) is again given by eq. 2.10. 

In the left panel of Fig. 3 we show the particle spectrum n(x, l,t,0) for impulsive 
injection for Kraichnan turbulence and for A = 0.3 and {^Bq/ L^pm) = 0.2. All timescales 
(rates) are in units of l/ko (ko). For early times {t < 0.1), the Green's function behaves 
similarly as in the cases considered before. For intermediate times (0.1 < t < 1), however, 
the acceleration rate has started decreasing while the adiabatic loss rate stays constant. 
For even later times, escape becomes dominant since its rate is cx 5^ cx exp [2Xt\ and the 
Green's function becomes heavily suppressed. This behaviour is even more prominent in the 
spectrum for steady injection, see right panel of Fig. 3. We show the particle spectrum for 
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Figure 4. Particle spectrum n{x,l,t,0) for impulsive (left panel) and steady injection (right panel) 
for Kraichnan turbulence. We assumed — ;> oo, a{t) — aoexp[— t/A] and k(t) = A:oexp[— ^/k] with 
flo = —0.5, X — 2, K = 1. The solid lines are for fixed times t = 0.01, 0.03, 0.1, 0.3, 1, 3 and the dashed 
line denotes the steady state spectrum. All timescales (rates) are in units of 1/fco {ko). 



intermediate and late times (t ^ 1) only since for early times (t < 1) it looks very similar to 
the previously considered examples. Around t ~ 5 the spectrum starts becoming noticeably 
asymmetric because of the dominance of adiabatic losses over acceleration. The escape rate 
l/xd oc g'^{t) oc exp [2At] keeps increasing and particles injected early have already escaped. 
However, particles which are injected late {t > 10) do not get accelerated anymore and hardly 
have time to lose energy before escaping. The resulting spectrum therefore converges against 
the 5-like injection though with decreasing amplitude. 

3.4 No escape 

Assuming a vanishing escape rate, Td — )• oo, the Ricatti equation can be directly integrated, 
= (g - 2fa^kg a = ao (^1 - {q - 2fao dt'kgj . (3.20) 

Furthermore, ^l^{t) and (p{t) are 

V'W=(^y and ^{t)={^{q-2f j\t'kg^ - . (3.21) 

For illustration, we consider both k{t) and a{t) to be exponentially declining, i.e. 
a{t) = aoexp[— t/A] and k{t) = fco exp [— t/zt]. In Fig. 4 we show the particle spectrum 
n(x,xo,t, to) for impulsive and steady injection for Kraichnan turbulence and choose oq = 
—0.5, A = 2, K = 1. All timescales (rates) are in units of 1/fco (fco)- At early times, the 
Green's function (Fig. 4, left panel) again behaves very similarly to the cases considered 
above. For late times, however, it stalls as both a{t) and k[t) approach zero. For momenta 
X > 1 and steady injection, the steady state is reached rather quickly as only particles in- 
jected early enough can reach higher momenta. Close to the injection momentum however, 
the spectrum still changes even at late times due to the injection of new particles. As t — ?■ oo, 
the spectrum reaches its steady state at all momenta except the injection momentum xq. 
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4 Summary and conclusion 

We have presented a new solution to the Fokker-Planck equation for stochastic particle 
acceleration by plasma wave turbulence. In extension to previously known solutions we 
allow the rates for stochastic acceleration, both diffusive and convective escape as well as 
adiabatic losses to be time-dependent. Furthermore, we do not need to constrain ourselves 
to specific values of the turbulence spectral index q and can therefore apply this result to 
the phenomenologically interesting cases of Kolmogorov {q = 5/3) and Kraichnan [q = 
3/2) turbulence. We have investigated four examples to illustrate the qualitatively different 
behaviour of spectra due to the time-dependent rates and extended range m q. In the first 
example we constrained ourselves, however, to constant rates and neglected adiabatic losses 
in order to compare our solution to previously presented results which constitutes a non- 
trivial test of our calculation. We have also considered the case of constant acceleration 
and exponentially decreasing adiabatic loss rate which leads to rather hard spectra. In 
contrast, with exponentially decreasing acceleration and constant adiabatic losses, the spectra 
become very soft and no steady-state solution exists. Finally, for infinite escape time we have 
presented an example with exponentially decreasing acceleration and adiabatic loss rates 
which leads to rather soft spectra in particular close to the injection energy. 

A few words about boundary conditions are in order. It was pointed out [41] that the 
initial value problem (IVP), eq. 2.1, with boundary conditions at x = and x — )• oo is 
singular (see Ref. [41] for a definition and discussion of a singular IVP). Therefore, the status 
of the boundary conditions is not clear but must be determined more carefully. Furthermore, 
the spectral theory of second order differential equations proves helpful in investigating the 
conditions under which a steady state solution exists. Unfortunately, it is not possible to 
analyse the problem at hand in this framework, since the variable coefficients of the partial 
differential equation prevent the formulation of an equivalent boundary- value problem (BVP) 
of Sturm-Liouville type. However, the transformed IVP, i.e. the heat equation type eq. 2.11, 
must satisfy boundary conditions for / at p — )• 0, oo similar to those that the original IVP, 
eq. 2.5, satisfies for / at x — >• 0, oo since both are connected by non-singular transformations. 
We find that /) — t- 0,oo are both limit points (for a definition see Ref. [41]) such that the 
appropriate boundary conditions at p — )■ 0, oo are simply ||/|| < oo which justifies the choice 
of the bounded solution in deriving eq. 2.11. Finally, we note that the discussion about 
the existence of the steady state solution cannot be applied here either due to the variable 
coefficients. As we have seen in the above examples, the existence of the steady state is largely 
determined by the time-dependence and asymptotic behaviour of the acceleration and loss 
rates, k{t), a{t) and T^{t) and therefore needs to be discussed on a case-by-case basis. 
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